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Abstract 

We show that longwavelength excitations of the quark-gluon plasma are de- 
scribed by simple kinetic equations which represent the exact equations of motion 
at leading order in g. Properties of the so-called "hard thermal loops", i.e. the 
dominant contributions to amplitudes with soft external lines, find in this approach 
a natural explanation. In particular, their generating functional appears here as the 
effective action describing long wavelength excitations of the plasma. 
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Significant progress has been achieved recently in understanding long wavelength 
excitations of a quark-gluon plasma|l], ^ [|. In equilibrum, at high temperature, such a 
plasma may be viewed as a gas of weakly interacting, massless quarks and gluons. When 
coupled to weak and slowly varying perturbations, this system may acquire a collective 
behaviour on a length scale ~ 1/gT, where T is the temperature and g the coupling 
constant, assumed to be small. In this letter, we present a consistent and physically 
intuitive description of such long wavelength phenomena, based on a set of coupled mean 
field and kinetic equations, to which the exact equations of motion reduce in leading order 
in<7g|. 

The kinetic equations encompass all the so-called "hard thermal loops" and clar- 
ify the nature of their remarkable properties, left largely unexplained by their original 
derivation in terms of Feynman diagrams^], [7], |2[. These hard thermal loops (HTL) are 
the dominant corrections, at high temperature, to amplitudes involving soft external lines. 
(Following the usual terminology, we call an energy or a momentum "soft" when it is of 
order gT, and "hard" when it is of order T; at equilibrum, most particles are hard.) As 
they are of the same order of magnitude as the corresponding non-vanishing tree level 
amplitudes, HTL need to be resummed consistently in higher order calculations Jl|, [2|]. 

Our equations isolate consistently the dominant terms in g in the hierarchy of equa- 
tions which describe the response of the plasma to weak and slowly varying disturbances, 
i.e. varying on a scale of order 1/gT. In doing so, we treat bosons and fermions on the 
same footing and introduce an average fermionic field ip(X) in parallel with the average 
gauge field A^X). A noteworthy feature of the present problem is that g, besides mea- 
suring the interaction strength, controls the wavelength of the soft space-time variations 
and, for consistency, also the strength of the mean fields. One finds for example that, 
in order for the deviations away from equilibrium to stay small, the gauge field strength 
tensor should be at most of order gT 2 ; then gA ~ gT is of the same order as the derivative 
of a "slowly varying" quantity [|], H|. 

The dominant interactions which determine the response of the plasma are those 
which take place between the hard particles and the soft mean fields. In leading order, we 
can neglect the direct interaction between the hard particles, which allows us to truncate 
the equations at the level of the 2-point functions. The motion of a given hard particle 
is only slightly perturbed by its interaction with a soft mean field. However, because the 
mean fields vary over distances much larger than the interparticle distance (~ ^-/T), they 
affect coherently many hard particles, giving rise to collective "polarization" phenomena. 
These show up as "induced sources" which add to the external ones in determining the 
properties of the mean fields. 
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The equations for the quark and gauge average fields are 

ip x if)(X)=r,(X)+r, ind (X), 



(1) 



and 

[rr x , F„„(x)] a - g^x^t^ix) = j;(x)+j; nda (x). (2) 

Here 7] and are external sources, a = 1, . . . , A^ 2 — 1 are color indices for the adjoint 
representation of the SU(N) gauge group, while //, v — 0, . . . , d — 1 are space-time indices 
(d — 4 throughout this work). The covariant derivative is _D M = <9 M + igt a A^ and i 7 ^ = 
[D^, D v ]/{ig). When using a covariant gauge, one should also consider equations for 
ghost mean fields. These are not written here as it turns out that they are trivial, i.e 
there are no induced sources for the ghost mean fields ||. After functional differentiation, 
the induced sources yield the one particle irreducible amplitudes with soft external lines. 
Thus for example the fermion self-energy is given by E(a;, y) = 5r] md (x)/5ip(y), the gluon 
polarization tensor by U°^ u (x,y) = 8j™ da (x)/5A'j! l (y), etc... 

In leading order, the induced sources can be expressed entirely in terms of 2-point 
functions. For example, r] ind (x) = g^ v t a {Af / {x)ifj{x)) c , where the subscript c indicates 
a connected expectation value which, in leading order, involves only the hard particles. 
The abnormal quark-gluon propagator which enters rj tnd is nonvanishing only in the pres- 
ence of the fermionic mean field if). It vanishes in equilibrum, as does the induced color 
current j . This current receives contributions from both fermionic and bosonic par- 
ticles, and accordingly may be written as j = j{ + jb- The quark contribution is 
jf a = g(^(x) r y lJ 't a 'il)(x)) c . In an arbitrary covariant gauge, the bosonic piece involves con- 
tributions from both hard gluons and hard ghosts ||. The polarization phenomena may 
be induced either by a gauge field A, or by a fermionic one, if). Correspondingly, we set 
jt = jf' + jf and jb = jh + 3b ■ Here, for instance, jf is the color current associated to the 
collective motion of hard fermions induced by the soft mean field A. It is independent of 
the fermionic field if). In contrast, we shall see that has generally a dependence on A 
imposed by gauge covariance. Note also that the ghosts do not contribute to j^, as they 
have no direct interaction with the fermionic mean field. 

In order to implement the condition that the average fields are slowly varying, it is 
convenient to use the Wigner transform of the 2-point functions, such as 

S(k,X) = I d ±se ik - s {Tif;(X+^(X-^) c , (3) 

for the quark propagator. Note that, in contrast to other authors, we do not insist on 
defining manifestly gauge covariant Wigner functions M. Covariance will be recovered 
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later, when calculating physical quantities. The induced sources are then integrals over k 
of Wigner functions, which we refer to as A;-space densities. For example: 

V md (X)=gJ-0^t a YIC:(k,X) J (4) 

where K,®{k,X) is the Wigner transform of the quark-gluon propagator K®(x,y) = 
(Tip(x)A a u (y)) c . Thus l£(k,X) = t a ljC a (k,X) = t a YK( k , x ) is the fc-space density for 
rf nd {X). The densities for the induced currents are denoted by J~ { A h (k, X) and J^ h {k, X). 
For example, the fermionic density is J^ik-, X) = r Ti r y ll t a S{k 1 X), where the trace refers 
to both spin and color indices. 

At leading order in the coupling g, the equations for the densities of the induced 
sources read|| 

(k ■ D x ) X) = -i g -(d - 2)C f {A{k) + &{k))W>{X), (5) 



k-D x ,J*{k,X) =2gN f k,kPF; u d»A(k), (6) 



k-D x , Jfl(k,X)] a = igk^{^X)t a Uk,X)-lfi{k,X)t a ^X)} 

- gk,f abc ^(X)t b lf c (k,X) -Tfi b (k,X)t c ^(X)} , (7) 

k- D x , J^(k,X)] a = gk^f^ {TP(X)t b )£ c (k,X)-lft b (k,X)t c iP(X)}, (8) 



k-D x ,J h \(k,X)\ =gN(d-2)k,k»F^crA(k). (9) 

In these equations, T/C(k,X) = ^(/c,X)7°, Cf = (N 2 — 1)/2N is the quark Casimir, 
Nf is the number of quark flavors and f abc denote the structure constants of SU(N). 
Furthermore, A(k) = p (k)N(k ) and A(k) = p (k)n(k ), where po(k) = 2ne{ko)5{k 2 ) is 
the spectral function for free massless particles, and N(k ) and n(k ) denote respectively 
boson and fermion occupation factors. The factor (d — 2) in Eqs. ([|) and ([]) reflects the 
fact that only the transverse gluons effectively contribute to the densities. 



Eqs. (g-|y) have a number of interesting properties: i) They are independent of the 
gauge fixing parameter A which enters calculations in general covariant gauges j|, |5|. ii) 
In their right hand sides, all possible vacuum contributions cancel, iii) They transform 
covariantly under a local gauge transformation of the mean fields A^, ip and ip. The 
densities jf , and If, involve Wigner transforms which are gauge covariant in leading 
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order. The current induced by a gauge field involves fc-space densities, J A and which 
are derived from non covariant Wigner functions. However, these densities are defined 
up to a total derivative with respect to k which does not contribute to the integrated 
current. We have used this freedom in order to make the densities explicitly covariant 
@, || . iv) The symmetry between Eqs. (|6]) and @ reflects the fact that hard quarks and 
gluons respond similarly to a soft gauge field. The same symmetry is apparent in Eq. @ 
expressing the effect of the fermionic mean field on the hard particles, v) Note finally the 
presence of the factor 5(k 2 ) in the r.h.s of these equations. This reflects the elementary 
dynamics of the hard particles: they remain on their unperturbed mass shell and only 
undergo essentially forward scattering on the mean fields. 

As these remarks strongly suggest, the motion of the hard particles described by 
Eqs. ([§-[5p exhibits many features of classical dynamics. This becomes more transparent 
if one makes explicit the structure of the various densities implied by these equations. For 
instance, Eq. @ implies 

J f \{k,X) = 2k fl N f t a {27r5{k 2 )) [e{k°)5n a + {k,X) + e{-k°)5n a _{-k,X)] , (10) 

where 5n± = 5n±t a are fluctuations in the quark color densities induced by the gauge 
field. These fluctuations satisfy (with = |fc|) 

■ D x , 8n ± (k, X)] a = T gv ■ E a (X)^-. (11) 

In the abelian case, this equation coincides with the linearized Vlasov equation. Here, 
the color electric field not only modifies the motion of the particle, but also induces a 
"precession" of the densities in color space. The other equations may be given similar 
interpretation. Thus, Eq. (j5|) for the quark-gluon Wigner function Iffi describes fluctuations 
where, under the action of a soft fermionic mean field, quarks are converted into gluons 
and vice-versa. 

The total current induced by a gauge field A is j A = jf + j£. Its fc-space density, 
J A = J t A + J A , satisfies 

~k ■ D x , J A (k, X)]=g k^k- F(X) ■ d (2N f A(k) + N(d - 2) A(k)) . (12) 

This equation generalizes the Vlasov equation to nonabelian plasmas. Previous attempts 
to derive such an equation led to more intricate results. However, they were based on 
different approximation schemes which mix leading and non leading contributions in g and, 
as such, are not entirely consistent^]. One can also combine Eqs. (0) and (^) into a single 
equation for the total current density induced by the fermionic fields, = jf + jt ■> 

k-D x , J*(k, X)} = ig k^ t a {$(X) t a X) - TfC(k, X) t a ^j(X)} . (13) 
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This equation is similar to the corresponding one in the abelian case M. In doing the 
sum of Eqs. ([?]) and (H), the typical non-abelian effects cancel; these are contained for 
example in the second braces in Eq. (|7|), and involve the 3-gluon vertex leading to gauge 
field insertions on the hard gluon lines. This kind of cancellation was first noted by Taylor 
and Wong || in relation with the HTL's for amplitudes involving one pair of quarks and 
any number of soft gluons (albeit their proof is only explicit up to three external gluons). 

We have thus reduced the set of equations (|§-|9|) to three fundamental equations, 
namely Eqs. (H), (0) and (p~3f) for the densities of the induced sources. From the gauge- 
covariant character of these equations, it follows that, under local gauge transformations, 
7] md transforms like ip, while j md transforms like F^ v . Provided the external sources are 
chosen so as to satisfy the same property, the mean fields equations (|l|) and @ are then 
gauge covariant, as are the classical equations of motion derived from the QCD action. 

Eqs. (|5|,|1^,|13|) contain all the information on the generalized polarizability of the 
plasma. After solving Eq. (§), we compute the fermionic induced source according to 
Eq. (H) and obtain 

v ind (X) = -ioj 2 \ —i \ duU(X,X -vu)ip(X -vu) 



An Jo 

= J d 4 Y6i: A (X,Y)ij(Y), (14) 

with = Cf(g 2 T 2 /8) and = (1, v), v = k/et- The angular integral runs over all 
directions of v, and U(x, y) is the parallel transporter along a straight line joining x and 
The kernel 5Y>a is the self energy of a soft fermion propagating in a background 
gauge field. The current induced by fermionic mean fields results from Eq. (|i~3|): 

j*(X) = gu 2 t a j —v,v u j o dtj Q ds 

ip(X - vt)YU(X - vt, X) t a U(X, X - vs)^j(X - vs) 
= gt a J ' d^Y^iY^Y^X-Y^Y^iY^. (15) 

The correction 5T A to the quark-gluon vertex may be easily read out from this equation. 
Finally, the current induced by soft gauge fields is determined from Eq. ([12]) to be 

—v. / du U(X, X - vu) FqAX - vu)v j U(X - vu, X). (16) 

4:71 JO 

Here ioi = (g 2 T 2 /9)(N + Nf/2) is the plasma frequency. By successive functional differ- 
entiation with respect to A of Eq. (|l^) we derive corrections to the equilibrum amplitudes 
for soft gluon fields. All the amplitudes obtained in this way, as well as those contained 
in 5T,a and in ST A , coincide with the HTL's of the diagrammatic approach. 
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Eqs. ( |I4} - |T6| ) imply the following covariant conservation laws for the induced currents |J: 



and 
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D»>tf =0, (17) 



igt a [ft a rf nd - fj ind t a 



By differentiating these equations with respect to the fields, one obtains relations between 
HTL refered as "QED-like Ward identities" in Ref.@, g. Finally, by using these relations, 
together with the Jacobi identity [_D M , [D u , F utl \] = 0, we see that the external sources 
must satisfy 



[D»,j„}=igt a (ipt a r}-f)t a TP). (19) 

in order for the mean fields equations (|T],0) to be consistent. 

By eliminating the induced sources from Eqs. (110), using their explicit expressions 
(|b|-[l6|), one obtains nonlinear equations of motion which generalize the Maxwell equations 
in a polarizable medium. In particular, for vanishing external sources, these equations 
describe the normal modes of the plasma. Note that in general, as a consequence of gauge 
covariance, quark and gluon modes mix. 

These non linear equations for the mean fields can be generated by the minimal 
action principle applied to an effective action S e ff = So + Si nd . Here So is the classical 
QCD action, while Si n d contains the effects of the interactions between the soft fields and 
the hard particles of the plasma. It follows that Si n d must satify 5Si n d/Sip(X) = rj md (X) 
and 6S ind /5A%(X) = jf da {X). These conditions are satisfied by S ind = S f + S g , with [g 

dtt 



r firs I r f — 
S f = -iul J — J <fX ^ dui/;(X)^U(X,X -vu)i/;(X - vu), (20) 



and 



Sn = l^f — ( d 4 X r du / du' 



Ju 



2 P J 4tt 

tr |>F M (X) U(X, X - vu') v u F uX (X - vu') U(X - vu', X)} , (21) 

where the trace acts on color indices only. This gauge invariant action coincides with the 
generating functional for HTL's derived in |J on the basis of gauge invariance. Here S e ff 
has a different, more physical, interpretation: it is the classical action describing long 
wavelength excitations in the hot quark-gluon plasma, at leading order in the coupling g. 
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